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Abstract 


We consider an inverse boundary value problem for the heat equation 
dtv = div x (yV^u) in (0, T) x P, where Pis a bounded domain of R 3 , the heat 
conductivity 7 (t, x) admits a surface of discontinuity which depends on time 
and without any spatial smoothness. 

The reconstruction and, implicitly, uniqueness of the moving inclusion, 
from the knowledge of the Dirichlet-to-Neumann operator, is realised by a 
dynamical probe method based on the construction of fundamental solutions 
of the elliptic operator —A + t 2 -, where r is a large real parameter, and a 
couple of inequalities relating data and integrals on the inclusion, which are 
similar to the elliptic case. That these solutions depend not only on the pole 
of the fundamental solution, but on the large parameter r also, allows the 
method to work in the very general situation. 
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1 Introduction 

1.1 Inverse heat conductivity problem 

Let P be a bounded domain in R , with lipschitzian boundary T = <9P, and consider 
the following initial boundary value problem 



divj, (”/X7 x v) in Pt = (0,T) x P, 
/ on r T = (0,T) x r, 

Vo on P, 


( 1 ) 


where 7 = € JF 1 ,oo ((0, T); L°°(P)) with the following properties: 

(C- 7 ) There exist a positive function (t,x) 1 —> k(t 7 x ) and, for all t £ [0,T], a 
non empty open set D(t) C P, such that 

• < k < C in Dt := U[ 0 .T]{i} x D{t) for some C > 1, 

• k — 1 has a constant sign in Dt , 
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~( t x) = { 1 if * ^ D{t) ' 

^ ’ ' \ k(t,x) if x £ D(t). 

We don’t assume any smoothness neither on D(t) nor that dD(t) D P = 0. 

Our main purpose is to study discontinuous perturbations, however, we allow 
7 (i, •) to be continuous. Hence, we impose the following assumption. 

(C-D) inf x ^k | k(t,x) — 11 > 0 for any compact set K C D(t), for all t £ [0,T]. 

We shall consider a large parameter r > 0 and allow the initial data vq(x) to 
depend on r, under the following condition: 

(C-0) There exist r-independent positive constants C, Iq such that 
||w||z, 2 (n) < Ce Tl °, for all r. 


Physically, the region D(t ) corresponds to some inclusion in the medium with 
heat conductivity different from that in the background domain 0. The problem we 
address in this paper is to determine D(t) by using the knowledge of the Dirichlet- 
to-Neumann map (D-N map) : 

A 7l „ 0 : / ^ d v v(t,x), (t,x)£T T , 

where v = V(y;/) denotes the unique solution of ( 1 ), v is the outer unit normal 
d 

to T, and d v = — = v • V x . In physical terms, / = f(t,x) is the temperature 
av 

distribution on the boundary and A 7i „ 0 (/) is the resulting heat flux through the 
boundary. 

The above inverse boundary value problem is related to nondestructive testing 
where one looks for anomalous materials inside a known material. 

To clarify our purpose, we remind briefly Ikehata’s probe method for the elliptic 
inverse problem. 


1.2 The elliptic situation 

In the probe method for the well-known elliptic situation, Problem (1) is replaced 

by 

J div x (yV^u) = 0 in U, , . 

\ v = f on T, [ > 


Dt is replaced by an open set D C fl. The Dirichlet-to-Neumann operator A 7 
is a mapping: H^(T) 5 f d v v £ where v is the unique solution 

of (2). The probe method (see [9]) starts by considering the fundamental solution 
ho{x) = 47r |^_y| of — A/io = 6 y , with pole y £ VL . Then, one approximates ho outside 

a needle E C fl with one end on T and the other one being y by a sequence {hj}j> \ 
such that — A x hj = 0 in f 2 , and estimates f D \Vhj(x)\ 2 dx (or f D \'S7h 0 {x)\ 2 dx) 
thanks to the following couple of inequalities: 


/ Whj(x)\ 2 < [ (A 7 (/ij| r ) — d v hj)hj\r 
Jd Jr 


— C \Vhj(x)\ 2 , 
Jd 


1 

C 


( 3 ) 
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for some C > 1 which does not depend on hj. However, in the parabolic situation, 
inequalities like (3) are unclear, except in the static case D(t) = D( 0), for a reduced 
class of functions as h(t,x ) = e T t p{x) : where r > 0 is a large parameter and p 
satisfies — A p + r 2 p = 0 in fh See [7], and [ 8 , 10] also for a similar approach. 

In our work we built fundamental solutions for parabolic operators which slightly 
differs from the heat operator C\=dt — A x or to its adjoint C\ = —dt — A x . Then, 
we obtain a couple of inequalities corresponding to (3), and, thanks to that, we run 
the probe method for the reconstruction of Dt- 

1.3 Special fundamental solutions 

Let t > 0, y £ R 3 , we consider the following function 

Pt, v (x) := -e~ Tr , r = \x-y\. 

It satisfies p T}V £ L 2 (R 3 ) and 

— A p T + r 2 p T = S(x — y) in R 3 . (4) 

Thus, the functions h T (t,x ) = e T t p T (t,x ) and h*(t,x) = e~ T t p T (t,x ) are respec¬ 
tively solutions to the following parabolic equations 

C 1 h:=9 t h-A x h = S{x-y), C\h* := -d t h* - A x h* = S(x - y). (5) 

Of course, since the right term in Equalities (5) is not 6(t— s)(gi8(x—y), our function 
h T differs from the usual fundamental solution G SiV (t , x) = e " 4(t ~ s) , with 

parameter (s,y) elx R 3 , considered in [1, 2, 4]. 

Let Q' be any smooth bounded domain containing O. Let E be an finite or 
infinite globally lipschitzian curve in D' with parameter y = y{t ), t £ M, such that 
y(t) VL for t < 0 . We shall use the restriction E|_i< t <T+i only, which plays 

the role of a Ikehata’s needle, but starts at point y(— 1) outside H and can be 

self-intersecting. We put 

Pr,y(t) =■ Pr, s(t, x) , (t, x) £ R X R 3 , 

and replace the functions h T ,h*, respectively by solutions of 

(£i + rr](t))U(t,x) = S(x — y(t)) in R x H, ( 6 ) 

(£* + TT](t))U*(t, x) = S(x — y(t)) in Kxd, (7) 

where functions 77 (f) are not necessarily the same in ( 6 ) and (7), but belong to 
L°°(R) with a r- independent upper bound: 

\v(t)\<P, t£R, (8) 

for some y > 0 that we shall precise later. (By q we denote the derivative of a 
function q according to the time variable t). 

The solutions U and U* we look for should be respectively written 
but we simply denote them by U T , U*. We construct them so that they have the 
following form: 

U T (t,x ) = e T t u T (t,x), XJ*{t,x) = e~ T t u* T (t 1 x) in f It, 
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with functions u T (t, x), u*(t,x ) sufficiently close to p Tt s(t,x) in the sense below. 

Lemma 1.1. There exists a function < p(t,X',r) G C([0, T]; (Mr)) satisfying the 

following points. 


;r) = l, feK, r > 0. 

There exists T\(p) such that for all r > n, we have 

< ip{t,X]r) < C(p) (t,x) elxfi, 


C(p) 


(9) 

( 10 ) 

( 11 ) 


| V x ip(t,x;r)\ < C(p), (i,i)elx(l, 
for some C(p) > 0 which does not depend on r or ( t,x ) G R x 12. The function 

U T {t,x) := e T t u T (t,x) = e T t ip(t,x-,T)p T ^(t,x) (12) 

satisfies (6) for all t > 0. 


Remark 1. Similarly we construct ip*(t,x;r) G C([0,T]; Hf oc (Mr)) such that the 
function 

U*(t,x) = e~ T t u*(t,x) = e~ T t ip*(t,x\T)p T ,T,(t,x) 
satisfies (7), with Relation 9 and Estimates (10), (11) where p is replaced by ip*. 

We choose functions p as follows. Let us consider the following function 

K(t) :=e- T ^- e \ teR, (13) 

where the parameters p > 0, 0 G [0,T], are independent of t and will be precised 
later. For U T , we put simply p = 0. For U*, we put 

p(t) = —— = /isgn (t - d). 

TK 

(Thus each p satisfies (8) obviously). 

1.4 Needle sequences 

Let T' G (0, T], we put Ht> = {(t,y(t)), 0 < t < T'}. We consider needle sequences 
{U* T ,}j , respectively associated to ( T,t,U t ,p = 0), (JCt',U*,p — ~Tk)- 
They satisfy: 

(jCi + rp(t))Uj = 0 in 12 t, 

Uj( 0,-) = U T ( 0,-) in 12, 

(Cl+rp(t))Ul T , = 0 in 12 T ,, (14) 

Ul T ,(T',.) = U* t (T',-) in 12, (15) 

and, for all open set V C 12 T ' such that d(V, S T /) > 0, 

(Uj,U* tT ,) (U T , U*) (strongly) in (H^(V)U H^ l (V)) 2 . (16) 
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Here, we use the notation 


H S ’P(V) = C^{Vy " 3 ' p , s,p> 0, 
with ||u|| s ,p = || \\u{t, •)llffP(R 3 ;dx))llff‘ i (R;dt)- 

Remark 2. Since y{ 0) ^ Q, thanks to (16), we have 


Uj,T' (0, •) ” f —>° C/*(0, •) (strongly) in L 2 (fl). 

1.5 Reflected waves 

Let Vj := V( 7 , Uj\r T ) be the solution of (1) with data / = Uj |r T , put Wj = Vj — Uj. 
By denoting £ 7 = d t — Va ;((7 — we have 


jC 7 Wj = V x ((j- l)V x Uj) in fl T , 
Wj = 0 on Tp, 

Wj\t=o = Vo -Ur in Cl. 

If St n Dt = 0 then, in view of (17), we have 

Wj j ^°W T in C([0,T],H 1 (n)), 

where W T is the unique solution of 

C^Wr = V x ((y — 1)V X U T ) in Sir, 
W T = 0 on Ft, 

W T | t=0 = v 0 -U T in Cl. 


Thus 


Vj 3 4° Vr = Ur + Wr ill H°’\V), 


(17) 

(18) 


(19) 


( 20 ) 


( 21 ) 


for all open set V such that V C fir \ £t- The function V T satisfies 

CjVr = S(x — y{t)) in CIt, 

V T = U T on Tp, 

K|*=o = iii ci. 

1.6 Pre-indicator sequence and indicator function 

Let 6 e (0, T), p > 0, and put 

Ij(r,n,6,T ) := f (d v Wj)U] T da(x)n{t)dt , 

JTt 

I 00 {t,P,6,T) := [ (7- l)V X V T V x U*r dxK(t)dt + [ [KWrU*r)l dx, 

where k is defined by (13), and da(x) is the usual measure on the boundary T. 

Remark 3. As opposite to dynamical probe methods in [2, 10], the indicator func¬ 
tion loo depends on the needle E. 
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The following result ensures us that we have the knowledge of 1^ (r, /X, 6 , T) from 
the Cauchy data of (1) when Ey does not touch Dt- 

Lemma 1.2. Assume that Et^iDt = 0. Then, for all 9 £ [0,T], y.,r > 0, we have 
Ij (r, /x, 0, T) ^°° /oo (r, ax, 0, T) £ K. (22) 

1.7 Main theorems 

The following theorem separates the cases E t D Dt = 0 and E t D Dt ^ 0- In the 
second case, we can put 

T* := min{T' £ [0,T]; E T , D D^ ^ 0} > 0. 


We put also 

A it', e ■= k max Q, > 0 <9 <T' <T, 

where Iq is the constant in (C-0). 

By d (y,X) we denote the ’’distance” from the point y to the set X, and by 
d(Y, X) the ’’distance” from the set Y to the set X: 

d(L,X) = inf d (y,X), d (y,X)= inf \x - y\. 

y£Y xGa 


Theorem 1. Assume (C-D) and let 6 £ (0,T). Assume that Et D Dt = 0. 

Then there exists > /XT.e such that, for all y > /xi, there exist C(y) > 1, T\ > y± 
such that, under assumption (C-0), we have 

— [ | 7 - 1| \S7 x p Tt x\ 2 dxKdt < \Ioo(t,ij,9,T)\ < (23) 

C( At) J fi r 

C(aO [ |7 - 1| \V x p T ,z\ 2 dxKdt, 

Theorem 2. Assume (C-D). Then the following points hold. 

• Let 9 £ (0,T). For all y > y±, under assumptions (C-0), E t f~l Dt = 0, we have 

lim — In | /oo(r,Ax, 0,T)\ = -2d(y(0), D(6)) < 0. (24) 

T —>OO 7" 

• under assumptions (C-0), E t D Dt = 0, we have 

limsup lim lim — In | /oo(r, At, 9, T)\ = — 2d(E-r, Dt) < 0. (25) 

#6(0,T) M - ^ 00 T 

• Assume that E t D Dt ^ 0. 

Then, under assumption (C-0), we have 

lim limsup lim lim — In | /^(t, y, 9, T')\ = 0. 

T'/'T* g e (o,T') r->°o r 


( 26 ) 
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So we can detect if St touches Dt or not. Consider the function 

F(T') = limsup lim lim lim — In | lj(r, p,, 9, T')\, 0 < T' < T, 

Q£.(q r p /^ Mt —>-oo j—yoo T 

defined in [0, T*[. The problem is that T* is defined from F itself. The following 
result gives a strategy to determine T*. 

Lemma 1.3. There exists S > 0 such that ifT' < T* then T' < T'+\F(T')\/5 < T*. 

In practice, if we have the a priori knowledge of 6, we can determine T* as the 
limit of the sequence 

(tn)nCN ■ to — 0, tn +1 — t n T |T(t n ) |/<5, XL ^ 0. 

If t n+ i > T for some n, this indicates that St does not touch Dt- We write 
T* = T + 0 in this case. 

Corollary 1.1. From the knowledge of A 7i „ 0 we can compute : 

(1) the maximal time T* of a given curve E. 

(2) the connected components of Dt \ Dt that touch { 0 } x T. 

In particular, the connected components of Dt \ Dt that touch {0} x T are 
completly characterized from A 7i „ 0 . 

Let us make some remarks on Theorem 1, Corollary 1.1. 

Remark 4. - We emphazise the fact that we don’t make assumption on D(t), 
0 <t<T. 

- From the knowledge of A 7) „ 0 , we get the information not only on d(Y,T ’, D(T')), 
when T’ < T*, but on d(y(Q),D(0)), 9 £ (0 ,T'), too. 


1.8 Literature review 

Many articles solve a version of the Calderon inverse problem for the heat equation. 
The biggest part of them assume that the unknown coefficient 7 , or the unknown 
inclusion D , do not depend on time t: see for example [8, 10 ] 

But they are very few results concerning the time dependent situation D = 
D (t). The authors in [3] proved uniqueness of Dt under the assumption that the 
inclusion D(t) is x-lipchitzian for all t. They used a proof by contradiction and is 
not constructive at all. A more recent paper for a similar question is [11]. 

A reconstruction method by a dynamical probe method is performed in [2] , but 
it works for the one dimensional spatial space only. In this case, an easier way is 
to put as input the trace of plane waves U T (t , x ) = e T t+TX at the spatial boundary, 
then to compute the solution V T (t,x) by using an ansatz and an energy estimate: 
see [5]. This approach works in the x-nrultidimensional case if D(t) is a convex set: 
see [ 6 ]. 

With the aim of working the dynamical probe method in the ^-multidimensional 
case with, the authors of [4] computed the reflected solution at points (t, x ) close 
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to the lateral boundary of Dt- But this way is painful enough to assume that 
D{t) = D does not depend on t. Moreover, since the data these authors consider 
are the traces on Tt of Runge approximations of the fundamental solution G SyV of 
the heat equation, this method requires some regularity of dD (more precisely, they 
assumed that dD is of class C 1,a , for some a £ (0,1)). 

In [1] the stability of the operator 7 1 — > A 7i o is quantified, under the assumption 
that Dt belongs to some class K, contained in ’l°°((0,T); W 2 ’°°(TL)) D W 1,oo ((0,T); 
W 1 ' 00 ^)). 

We think that our approach, based on special fundamental solutions and in¬ 
equalities (23), would be able to extend the frameworks of many of these articles. 

1.9 Outline of the paper 

The paper is organized as follows. Section 2 is devoted to the proofs and is devided 
in several subsections. In subsection 2.1, we give basic estimates for integrals on 
Dt with weight functions as (p T ^K,)(t,x). In subsection 2.2, we build the functions 
U T , U* by proving Lemma 1.1 of section 1.4. In subsection 2.7 we prove Lemma 
1.3. In subsection 2.5 we prove Theorem 1. In subsection 2.6 we prove Theorem 2. 
In subsection 2.8, we build the needle sequences Uj, Uj y T'- 


2 Proofs 

Notations 

We consider the usual Sobolev spaces H m ( Id), where V C R fc is an open set, and 
also -ff m ((Ti, X 2 ); H m '(V)) = H rn (T 1 ,T 2 ;H m (V)), for m,m! £ Z. Here, H° is the 
well-know space L 2 . 

The formal parabolic operator is £ 7 = dt — V x • (yV^-). 


2.1 Basic Estimates 

In this part we establish some basic estimates on W T etc. 

Lemma 2.1. a) Let Ocl 3 be an non empty open set and r > 0. Then we have 
f \Pr y y{x)\ 2 dx<C—^—e- 2Td ^°\ VyeK 3 (27) 

Jo T d{y, O) 


j \VxPr,y{x)\ 2 dx < Ct (l + * ) Vy £ O. 


(28) 


b) Let O C R 3 is a non empty open. Let y (jL O and d > d(y,0). Then there exists 
C(0,y,d) > 0, such that 


f \\7 x p T ,y{.x)\ 2 dx > C(G,y, d)r 2 e 2rd , Vr > 0. 

Jo 

Proof, a) We check easily: 

Io Mx)fdX = i(^ e ~ 2Trdx - Li l , 1 o^ TrdT 

e~ 2s ds = (2rd(y,C>))- 1 e- 2rdfe ’ 0) . 


( 29 ) 


< 1 / 

T J s>rd(t/,0) 
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Hence (27). Similarly we have 

/ |V xPt, v ( x )\ 2 dx = 

Jo Jc 


(rr + l) 2 


2 Tr, 


< 


O 

h 


e “' ‘ dx < 
2 


(rr + l) 2 n Tr o 


(4?r ) 2 r 4 ’ Jr>d(y,0) 4 ^ 4 


dr 


(s + 1) 2 - 2s 


«>rd(j/iO) S 


ds 


Ct ( 1 + _ ) e~ 2Td{ - v '°^ 


rd(y,0) / 

b) Let us prove (29). Observe that the open set Od := {x € O; |x — j/| < d} is non 
empty. Hence we get 

J(y) := / |Va;p T y (x)| 2 dx > / ^vr e-dx^r'e- 

do J Od ( 47r ) 1 

with c = J 0d (i^pdx > 0. 

Remark 5. It is clear that if y £ O, or if dO is lipschitz and y £ dO, then 


/ |V x p T>y (x)| 2 dx = +oo. 

Jo 


Lemma 2.2. Assume that Tt O Dt = 0. Let 8 £ [0,T]. Let d > d(y(0), D(6)). 
Then, there exists y\ > 0, ri > 0 such that, for all y > y\, t > t±, we have 

1 


[ -l\\V x p T ^(x)\ 2 dxK(t)dt <—(l 

Jn T T 


) e -2 rlHy{0),D(0)) t ( 30 ) 


Td(y(e),D(8)) 

[ |7 — 1| |V x p r> E(x)| 2 dx«;(t)dt > C—e~ 2Td 
J Qt 1 M 


(31) 


/or some C > 0 which does not depend on y or r. 
Proof. 

a) Thanks to (28) we have 


I ■= I |7 — 1 | \TIxPt,y.{x)\ 2 dxn{t)dt 

hi T 

< Ct J (l4-I- N l p -2rd(y(t),D(t)) 


e K(t)dt. 


to v rd(y(t), £>(*)), 

By assumption on Dt and E, we have 

d(y(8), D{9)) + S\t -0\> d(y(t), D(t)) > d (y(0),D(0)) -6\t-0\, te [0, T], 

for some S > 0 depending on Tt and Dt only. Fix yi > 26 and let y > yi- Put 
= dfa(flkpW) w h 


I < 


C . re -2rd(y(S),D(e))^ 1 + 

_|_2 1 g-2' r d(i/(S),D(6i)) e re( 2 ( 5-/j) ^ _j_ 


1 


r(d(y(0), D(0)) - 5e) J J\ t _ B \ <e 
1 


e r(2S-c\t-e\ dt 


< e ~ 2 Td(y( 9 ),D( 9 )) f g ^ + 


7"(d(Er, Dt) 
2 


dt 


l At -2d v Td(y(0),D(9)) 


) +0(e- 5 '^)}. 
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for some S' > 0. Estimate (30) is proved, 
b) Put 

J(t) := [ | 7 (t,x) - 1| \S7 x p T , v (t){x)\ 2 dx. 
J n 

Thanks to assumption (C-D) we have 


J(0)>S / |V*p TlI ,( fl )(a;)| 2 da;, 
Jn 


with S = liminf^gQ | 7 ( 0 , x) — 1| > 0. Then, thanks to (29), we have 


J(6) > 6'T 2 e~ 2Td , 

for some S' > 0. Obviously, </(•) depends continuously on t, since 7 £ L°°(Qt)- 
Hence if |t — 0\ < e and e > 0 is sufficiently small, we have 

J(t) > U'T 2 e~ 2Td . 

We then have 

I := f J(t)n{t)dt > f J(t)n(t)dt > CT 2 e~ 2Td f n(t)dt 

J Qt J\t—0 \<e J\t—6 \<e 

> C-e~ 2Td {l - e~ Ty,e ). 

/i 

for some C > 0. Taking /ri > 7 , this proves (31). 


Lemma 2.3. Let 9 £ (0,T). Then there exist S > 0, /zi > 0 such that for all 
H > H 1 there exist C(fi) > 0, n > 0 such that 

f n(0)\U\ 2 dx < C(ia)e~ 2Td{yW ’ DW)+5 \ t £ [0,T], r > n, 

Jn 

U = U T (0), t/*(0); (32) 

[ k(T)\u*(T)\ 2 dx < c(n) e - 2T(dw) ’ Dm+s \ r > 1; (33) 

Jn 

f K(0)\W T m 2 dx < C(n)e~ 2Td(v{t) ' D(t))+s \ t £ [0,T], n > /xi, r > 1; (34) 

Jn 

[ k(0)\W t ( 0)U* (0)\dx < C'(/i)e _2T ,-»(*))+«) s t £ [0,T], r > n. (35) 
Jn 

Proof, a) Thanks to Lemma 1.1, (10), to Remark 1, to (27), and reminding that 
y(0) ^ ft, we have 

f \U\ 2 dx < C f \p TM o)\ 2 dx < c T - l e- 2Td( - y{ - 0) ’ n) < C'. 

Jn Jn 

Since k(0) = e~ Ty6 , we then have 

[ n{f))\U\ 2 dx < C'e~ Ty9 . 

Jn 


(36) 
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Choosing pi > sup( t]X ) e n T 2 ^ v ff ' >,x ) , we obtain (32). 

b) Similarly, choosing pi > sup( t ,x)en T 2d (r-e) z) > we obtain (33). 

c) Let us remind that W T ((), x) = vq(x) — C/ T (0, x). Thanks to Assumption (C-0) 
and to (10), (12), (32), we have 

[ K (0)\Wr(0)\ 2 dx < Ce~ Tfl0 {e loT + T- 1 e~ 2Td{v{0) ’ n) ) < Ce - 0 ^-^- 9 ^. 

Jn 

Choosing pi > / i T ,e + supr t>x ' )& Q T 2 ^v<f)^ x ) we obtain (34). 

d) Estimates (32), (34) and standart inequalities imply (35). 


2.2 Construction of the special functions 


Let R > 0 such that C {x \ \x\ < R} = Br. Put z = x — y(t), u T {t , x) = f(t , z). 
We have 

d t u T = d t f + zV z f = d t f - ijV z f, 

{d t - A z )f - y(t)W z f + (t 2 + TT](t))f = S 0 {z). 

Let /(£, fc) = (27 t) - 3 / 2 / zgR3 fit, z)e~ xkz dz be the Fourier transform of fit, •). It 
satisfies 

dtf + (fc 2 - iky(t) + t 2 + Ty(t))f = 1. 


Put p{t) = f 0 r/(s)ds. We then have 

rt 


fit , fc) = 

(2 tt)- 3 / 2 /", 

J— C 

and so 


f(t,z) 

= (2tt)- : 


1 


87T 3 / 2 

it r (£, x) 

1 

87T 3 / 2 


(37) 


fit, z) = (2tt)- 3/2 [ ([ e^^-^+^+yW-vi^dk] e r 2 (*-t)+T(p(«)-pW)d S) 

J —co k£M. 3 / 


e ~k 2 s+ik(z+y(t)-y(t-s)) dfc ) e -r 2 s+r(p(t-s)-p(t)) 
/o V-lfceR 3 


z -T 2 s+r(p(t-s)-p(t)) e ^^l^- s -3/2 ds _ 


(38) 


/o 


The integral in (38) is convergent if and only if x ^ y(£). We set 

u T 

ip := -. 

Pr,S 


Let us prove (10). We put 


9ris,t,x) 
fc T (s, t, x) 


0 T{p(t-s)-p(t))+-^(y(t-s)-y{t))(2x-y{t-s)-y{t)) 


1 




8tt 3 / 2 


o-3/2 


and write 


m t (£, x) 



k T (s, £, x)g T (s, £, x)ds. 











Let us observe that 
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Pr,y(t)( x ) = / k T (s,t,x)ds, l£l 3 . 

Jo 

In fact, the right-hand side in the above relation belongs to L 2 (R 3 ), with a Fourier 
transform corresponding to (37) where g is replaced by 1. Moreover, it satisfies (4), 
since it is a time-independent solution of (6). Hence this right-hand side is p T , y (t)( x )- 

Since p(-) is lipschitzian, we can fix K € R such that r > \K\ > ||p||oo- Let us 
remind that p = rj and rj satisfies (8). Thus the above inequalities impose p < t, 
and, in fact, 

p « t, as r —> oo. 

Then, if K > 0, we can write, for x £ ft C Br, x ^ y(t), 

k T (s,t,x)g T (s,t,x) = k T _x(s, t, x)e 2rKs+K s g T {s 1 t, x)ds 
< C(p, K, R)k T _K(s,t,x). 


Hence 


u T < C(p,K,R)p T> E in fly. 

For K < 0 we obtain 

U T (t,x) > Pr-K,y{t){ x ) > C"(M> K, R)p T ,y(t) {x), 

where C'(p,K, fi) > 0. Thus (10) is proved. Let us prove (11). Let x G Cl C Br, 
x 7^ y(t)> We use the following changes of variable 

T T 

s = —e a , s := —— e~ a , b = cosha (a > 0), a = rrb ', b' = 6—1. 


2t 

We then have 


2 r 


1 2 2 

= ^372 J e~ r2s ~^s~ 3/2 g T {s,t,x)ds 

V2 J e ~ TrC ° Shae ~ a/2 9 r(s, t, x)da 


87r 3 / 2 
1 ,2 r 


' —OO 

/»oo 


8tt 3 / 2 

1 ,2t 


2r r 00 

( — ) 1/2 / e _Tr cosha {e~ a ^ 2 g T (s, t, x) + e“/V(3, t, x))da 
r Jo 

o t r 00 / 

( —) 1/2 y e _rrb ^cosh(a/2)(c/ T (,s, t, x) + g T (s,t,x)) - 


87r 3 / 2 

sinh(o/2)(g T (s, t, x) - g T (s,t,x))jda 


u + — u 
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with 

u+{t,x) 


1 , 2r, 

87T 3 / 2 

1 ,2r 


2r Z" 00 

(—) 1//2 J [g T (s, t, x ) + g T (s, t, x))e~ Trb cosh(a/2)da 
__(^)i/2 e — i, s) + flT (S, i, x))e— b ' C ° s ^ ) 2) d6' 

^ 1 - rr r(9r(s,t,x)+g T (s,t,x)) C ° sHa/2) -- a da 

Jo 


sinh(a) 


x r (t,a;) = 


87t 3 / 2 r 
Pr,y{t) t P + {ti x )i 

i o T r 00 

^(-) 1/2 J (ffr (s, t, x) - g r (s, t, x))e~ Trb sinh(a/2)da 

1/21 -r M s,t^- 3 As,t,^ ia/2) - da 


Stt 3 / 2 7- 


sinh(a) 


= PT,y(t)<P (*,®)- 


Estimate (11) in Lemma 1.1 is a consequence of the relation ip = ip + tp + and the 
folllowing Lemma. 


Lemma 2.4. We have, for r < R and r sufficiently large (t > C(p,R)), 

\ip~(t,x)\ < C(n,R)/y/f, 

\<p + (t,x)\ < C(n,R), 

< C(jjl,R). 

Proof. We observe that 0<s<l<sfora>0, and 
cosh(a/2) = y/(V + 2)/2 = J + 1, 
sinh(a/2) = \JVJ2 = 

t(s — s) = r sinh(a) = r\/b 2 — 1 = r^b'ib' + 2) = —=\fa* — + 2r, 

V t 


(39) 

(40) 

(41) 


a 2 r 


< - + 
T 


a, 


r(s + s) = r cosh(a) = rb = r + rb' = r H- . 

T 

Hence, for 0 < r < R, we have 

9 r{s,t,x) = exp| — r(p(t — s) — p(t)) + -^(y(t — s) — y(t))(2r + y(t) — y(t — s))| 

< exp^rs|p|oo + \y\oo(r + s|y|oo)) < e c(r+rs) < Ce 2Cr+c - < Ce c ^, 
where C = C(p,R). The same estimate for g T (s,t,x) holds: 

max(< 7 r (s, t, x), g T (5, t, x)) < Ce c ~. 


















We thus have 


\g T (s,t,x) - g T (s,t,x)\ < (r\p{t - s) - p(t - s)| + - s) - y(t)) 

•(2x - y(t - s) - y(t)) - ^={y{t - s) - y(t))(2x - y(t - s) - 


•ma x(g T (s,t,x),g T (s,t,x)) 

< (|p|ooT(s - s) + i|y|oo(2 r + |y|oos) + ^|y|L( s + «)) e 


c# 


a 2 r 


< C(t(s — s) + r + s + s)e ct <C( -b \ — \[a + r H—(rb—))e 


oC§ 


,a 


< C {—b r)e c ~, C = C(n,R). 


Consequently, since 


<p (t,x) = 


1 


V2tt Jo 
^ /*°° 


, , . ..sinh(a/2) da 

(g T (s,t,x) - g T (s,t,x)) . , , J e 




sinh(a) 


rr 

da 


we then have, for r > C = C(p, R), 

r OO 

| <p~(t,x)\ <C (—I -r)e^~^ c 

Jo T 


2 cosh(a/2) 


da 


rr 




1 


< cl ^ + JL )e ^-^ da 


< c'(- + J-). 


Thus (39) is proved. Moreover, since 


V x g T (s,t,x ) = ——y—— g T (s,t,x), 

2s 

\V x g T {s,t,x)\ < ^|y|oo|ffr(s,t,x)] < Ce 


Ch 


we then have 


/•OO 

\X/ r y~{t,x)\ < C \g T (s,t,x)+g T (s,t,x)\ 

Jo 


1 


da 


2 cosh(a/2) 


C 




1 


da 


C / < C". 


V^ + iv 7 ^ 7 

da 

\[a 


Thus the estimate of V r ip (t,x) in (41) is proved. Let us prove (40). Since 
cosh(a/2) 1 1 


sinh(a) 2sinh(a/2) \/2a’ 
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we then have 


1 r°° 

(p + (t,x) = — 7 =f (g T {s,t,x) + g T (s,t,x)) 

V 2.71 Jo 


cosh(a/2) _ a da 


sinh(a) 


rr 


and so 


\<p + (t 1 x)\ < C | / e c ~~ a _ da < C'. 

Jo V2a 


/ o V%a 

This proves (40). Similarly, thanks to (42), we have 

1 


\S7 x ip + {t,x)\ = 


C 


I 1 

/ (^xgr(s,t,x) + \/ x g T (s,t,x))—=e a da 
Jo \J2a 


OO 

e C '-“ 1 


.— da < C ', 
y/2a ~ 


which proves (41) for V x <£> + . So Lemma 2.4 is proved. 

By observing that = ip~ + (p + , Lemma 1.1 is then proved. 

2.3 Properties of reflected waves 

Putting 

*r = ( 7 -l )V X V T + V X W T =- 1 V X V T -V X U T7 
vj = (7 - 1 yVxVj + V x Wj = -fV x Vj - V x Uj, 

we have 

div x 4/ T = div x ('y(t,-)V x V T ) - A X U T = d t W T in n T , 
div^Tj = div x ('y(t,-)V x Vj) - A x Uj = d t Wj in Cl T . 

If St D Dt = 0 then, in view of (43), (16), (19), we have 

in C([0,T],(L 2 (L1)) 3 ). 


(43) 


(44) 

(45) 


2.4 Proof of Lemma 1.2 


Let us set T' =T which does not restrict the proof. We denote U* instead of U* T . 
Put 


Jj(t) ■■= [ (d v Wj)(t, x) U* (t, x)da(x), 

Jon 

Fj (t) := /( 7 -l )V x V 3 V x U*dx. 

Jn 

Thanks to (45) and integration by parts, we have 

Jj(t) = [ (v ■ ^j)U*da(x) = I T jV x U*dx+ f (di v x ^j)U^dx 

Jan Jn Jn 

= F 3 (t) + f (V X W : V X U* + d t W 3 U*)dx. 

Jn 
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Thanks to (18), (14), we compute 

/ S7 x Wj'V x U*dx = - I WjA x U*dx = f W 3 {d t U* + C\U*)dx 

= [ Wj(d t U* + —U*)dx = [ WjK-'dtinUfidx. 
J n K Jn 


Thus 


Jj{t ) = Fj(t) + f k 1 d t (WjU*K)dx. 
J n 


(46) 


Thanks to (14), (46), an integration by parts according to t brings 


Ii = 


Jj(t)n(t)dt= f Fj[t)n{t)dt + j [WjU^K^dx, (47) 

Jo Jn 


Thanks to (15), (19), (21), (47), and to Remark 2, if E t fl Dt = 0, then we 
have: 


*J 3^+oo 


[WjUj^Yo dx 


[w T u;^dx, 


[ F 0 {t)n{t)dt I T = f (7 - l)V x ^V x E7K(*)dt, 

J 0 J £It 

Ij ^4°° / oo (r,Ai,0,T) = j T + /[W^Klg’d®, 

J o 

which proves (22) in (1) of Lemma 1.2. 

2.5 Proof of Theorem 1 

For the sake of simplicity, we assume that 7 > 1. The case 7 < 1 is similar. We put 


Fr ■■= [ (7-1 )V x V T V x U*dx. 
Jn 


2.5.1 Lower Bound for 1 ^ 

Let us give another expression of 


I T = F T (t,x)n(t)dt. 
Jo 


We remind that 


T T* -2r frr 

U T = pe U T , p := 




where <p, p* are characterized in Lemma 1.1, and we put 

7 - 1 t, / ^iPlrrt 


ri := 


r 2 := 


-tt T (-^)E£, 

7 p 


7-1 


( 48 ) 


{(V^V + x p)p Tl Y,y xPt.y, + V x <£*V x <£|p Ti e | 2 j. 
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We firstly have 

( 7-1 )v x v T v x u: ( => l^la> T v x u; + l^lv x UrV x u; 

7 7 

= ri + Tzl^ T pV x U T e~ 2T2t + 1^±V X U T V X U; 

7 7 

= ri + I^ / 9 e - 2 r 2 t _ ^ r ( V;cWT ) / 3 e - 2 r 2 t + 

7 7 

= ri + —'I ,2 pe~ 2r2t - 4' T (V a; W T )pW 2r2 * + -—- | V x p T , s \ 2 tp*ip + r 2 . 
7 7 

Explanations - (1): write = ^(tE 1 + V X U T ). (2): use (48). (3): write (7 — 

1)V x 17 t = \I/ T — 7 V X W T in the second term. 

By integration by parts, thanks to (20) and (44), we have 


We thus have 

F t = 


Writing 


- / 'S’tV x W T pdx = / d t W T W T pdx + r 3 , 

Jn Jn 

r 3 := / 'b T W T V x pdx. 

Jn 


I — — -\V x p T ^\ 2 p*pe 2r t dx+ [ -ty 2 pe 2r t dx 
Jn 7 Jn 7 

7-1 


+ / (0tW T )W T pe dx + (-(ri + r 2 ) + r 3 )dx. 

Jn Jn 7 

f (d t W T )W T pe~ 2Tt ndtdx = f (t 2 p — ^-p)\W T \ 2 e~ 2Tt ndtdx 
J J Cl'j 1 

+ Uj v ' 


r2 e~ 2Tt K\n 


we then have 


Ioo(t,p,6,T) = I t + J n [W T U*K]^dx 
= fn T ^NxPr^^VKdtdx + J Qt 2 T pe~ 2T2t Kdtdx 
+ fn T (T 2 p— ip)|W T | 2 e^ 2T * ndtdx + \ f n (W 2 e~ 2rt )k | t =rdx + R, 
where we put 


R :— i?i + R 2 + i?3 + Ra + 7?5, 

Rj := / Tjndtdx , j = 1,2,3, 

J 



R5 ■= ~ f {W 2 K)\ t=0 dx, 

J o 


(49) 


(50) 


(51) 


We prove in the next Lemma that R is neglictable. 
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Lemma 2.5. There exist 5 > 0, p i > 0, so that for all p > p\, there exists t± > 0 
so that for all r > t\, the following estimates hold. 

|i?i| < r _1 f ^ 2 pe~ 2Tt ndtdx + CT~ 1 f \"f — 1\ \\7 x p T ^\ 2 ndtdx; (52) 

J J Qx 


I-R 2 I < Ct 1 I 7 -1| \V x p T ^\ 2 Kdtdx; (53) 

J Qt 

| 1 ? 3 j < f — T 2 pe~ 2r t ndtdx + C f \W T \ 2 e~ 2r f ndtdx] (54) 

JO.T 2 7 Jq t 

\R 4 \ <7 [ \W T \ 2 \ t = T e- 2TT K(T)dx + Ce- 2T(vW ' Dit))+5 \ (55) 

4 J n 

|i? 5 | < Ce~ 2T{v{t) ’ DW)+5) , t e [0,T]. (56) 


Proof, a) Since p T < r\\/ x p T ^\ and thanks to Remark 1, we obtain (52) 
and (53). 

b) Estimates (54) comes from (10), (11), from Remark 1 and standart inequalities. 

c) Thanks to Lemma 1.1, (10), to Remark 1, to (27), we have 

|i ? 4 |<7 f \W T \ 2 \ t=T e- 2TT K(T))dx + C [ \p T , v(T) \ 2 n{T)dx + \ \W r U^\ t=O K(0)dx 
4 Jn J n Jn 

d) Estimate (56) is (34). Then (55) comes from (32) , (35). This proves the Lemma. 
From Lemma 2.5, (51), (10), Remark 1, and since 

r 2 p- \p + C> C't 2 , C">0, 

for r sufficiently large, we obtain 

Ioo(t,p,6,T) > -^{ / -- -\X7 x p T ^\ 2 K,dtdx + [ ^ 2 pe~ 2T t Kdtdx 

^ 1 J q t 7 J n T 

+ f T 2 \Wr\ 2 e~ 2Tt2 ndtdx + f |W T | 2 | t=T e“ 2rT2 K(T)da:), 

for some C > 1, which implies the lower bound for \Ioo\ in (23). 

2.5.2 Upper Bound for I Q 0 

In the same way, putting 

r 6 := (7- ljV^V^^e- 2 ' 2 *, 


we have, 

( 7 - 1 )v x v T v x u; ( = } (7 - i)X7 x w t \7 x u * + (7 - 1 )v x u T v x u; 

= (7 - l)pV x W T V x U T e- 2T2t + (7 - 1 )V x U T V x Uf + r 6 
= 1 P \ 7 x W T \/ x U T e~ 2T2t - pS 7 x W T S 7 x U T e- 2T2t + (7 - 1 )V x U T V x Uf + r 6 
= ^pV x \W T \ 2 e~ 2T2t + P V x W T <b T e- 2T2t + (7 - l)V x U T V x Uf + r 6 . 
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Explanations - (1): write V T = U T + W T . (2): use (48). (3): write U T = V T — W T in 
the firs two terms. 

Thanks to (49), we have 

F T = — [ 'yp\/ x \W T \ 2 e~ 2r *dx - f d t W T W T pdx - r 3 

Jn Jn 

+ / ( 7 - l)V x U T V x U*dx + / r 6 dx. 

J£l J Q, 

Thanks to (49), (50), we obtain 

I T = — f 7 p\\7 X W T \ 2 e~ 2r t ndtdx — f (r 2 p — -p)\W T \ 2 e~ 2Tt ndtdx 

J Qj'jp J " 

[ (W?e~ 2Tt n) t= Tdx + [ (7 - 1)V X U T V x U*ndtdx 

J Jn Jn T 

—R 3 — 1?5 + Rq, 

where we put 

Rq := / rQn(t)dtdx. 

J Q'j-' 

Similarly to the above section §2.5.1, for r sufficiently large, we obtain 

Ioo{t,p,, d,T) < -^( [ (7 - 1)\\/ x p T ^\ 2 ndtdx - [ | V x W T \ 2 e~ 2T2t ndtdx 

7 t Jn T Jn T 

— I T 2 \W T \ 2 e~ 2Tt ndtdx — f | VE r | 2 | t=T e _2TT n(T)dx j, 

Jn T Jn J 

for some C > 1, which implies the upper bound for j/ool in (23). 

Thus Theorem 1 is proved. 


2.6 Proof of Theorem 2 

a) Assume that T < T*. Let 9 G (0,T) and p > pi, where p± is defined in Theorem 
1. Thanks to (31) and to Theorem 1, we obtain (24). Since in this formula the right- 
hand side does not depend on p, we let p tend to infinity. We then obtain (25). 

b) Assume that E t D Dt ^ 0, that is, T > T*. Let 0 < T' < T*. Thus, (25) 
holds with T replaced by T' . Taking the limit as T' /* T*, and since d(E 
tends to d(E-r*,-DT*) = 0 when, we obtain (26). 

2.7 Proof of Lemma 1.3 

First remind that F is defined at least for T' sufficiently close to 0. Observe that 
the function G: 

G : [0,T] D T' -2d(Y, T >,D T >) 

is lipschitzian, non-decreasing, and G(0) < 0. Thus, there exists 5 > 0 such that 

G{T”) < G(T') + 2<5(T" - T'), 0 <T' <T" < T. 

Hence, if G(T') < 0, then G(T r + \G(T')\/5) < 0. The result is proved since , G = F 
on [0, T*[. 
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2.8 Existence of Needle sequences 

Denote by £ c the open set ((—1, T+l) x iy)\£[_ 1T+1 ], with £[(,*'] = {(s, y(s)); t < 
s < t'}. Let us prove the existence of {uj}j. Let G(t, x) be the solution in T 2 ((0, T + 

l);fl o 1 (n / ))nJT 1 ((0 ) T+l);L 2 (n / )) c <7(1-1, T+l]; T 2 (Q')) of 

(C\ + Trj)G = 0 in Cl' T+1 , 

G|t=o = U T (P,-) in Cl', 

and d t G G T 2 ((0, T+l); iJ _1 (D')). Since U T € C([-l, T + 1]; L 2 (Cl’)) and d t U T G 
L°°((—1, T+l); TT —1 (ST')), then the function H = U T — G extended by 0 for t < 0 
belongs to C([— 1, T+l]; L 2 (Cl')) with d t H G T 2 ((—1, T+1); and satisfies 

(£i + rr])H (t, x) = X(o,+oo )(t) ® 5(x - y(t)) in (-1,T+ 1) x Cl', 

H\ t <o = 0 in Cl'. 

Let U be a simply connected open set in R 4 with lipschitz boundary, and such 
that U C Ti C . Observe that £ c \ U is connected. For all open set V such that 
U C V C V C Ti C , we have H\y G iL 1,2 (’F) , H\v(t,-) = 0 for — 1 < t < 0, 
and (Ci + rr])H = 0 in V. Similarly to [2, Theorem 2.3], there exists a sequence 
Hj G (— 1,T + 1) x Cl' such that 

(C 1 + Tri)H j (t,x) = 0 in (—1, T + 1) x Cl', (57) 

Hj\ t <o = 0 in Cl', 

and Hj converges to H in L 2 (U) (with the strong norm). We then have the existence 
of another sequence Hj G H^(T, C ) satisfying the homogeneous relation (57) and 
such that Hj converges to H in H loc (£ c ), as noted in [2, §3]. Finally, the sequence 
Uj = Hj + G satisfies Equation (57) in Cl' T , with Uj\t= o = U T (0) in O', and Uj -+ U T 
in H ioc i^T \ O H^l(Cl T \ St)- 

Remark 6. We used the fact that U T (0) G L 2 (Cl') only. Hence, the proof for 

the existence of U* in the dual case, where the triplet (C\,t,t = 0) is replaced by 

(C\, T — t,t = T), works similarly. 
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